Entanglement entropy of 2D conformal quantum critical points: 
hearing the shape of a quantum drum 
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The entanglement entropy of a pure quantum state of a bipartite system A U B is defined as the 
von Neumann entropy of the reduced density matrix obtained by tracing over one of the two parts. 
Critical ground states of local Hamiltonians in one dimension have entanglement that diverges loga- 
rithmically in the subsystem size, with a universal coefficient that for conformally invariant critical 
points is related to the central charge of the conformal field theory. We find the entanglement 
entropy for a standard class of z — 2 quantum critical points in two spatial dimensions with scale 
invariant ground state wave functions: in addition to a nonuniversal "area law" contribution pro- 
portional to the size of the AB boundary, there is generically a universal logarithmically divergent 
correction. This logarithmic term is completely determined by the geometry of the partition into 
subsystems and the central charge of the field theory that describes the equal-time correlations of 
the critical wavefunction. 

PACS numbers: 03.67.Mn, 11.25.Hf 



A major challenge in the theory of quantum critical 
phenomena is to understand those aspects that do not 
appear in the classical theory. The entanglement entropy 
at criticality is an example: a pure quantum state of a 
bipartite system AU B can become a mixed state, with 
an associated entropy, when restricted to subsystem A 
or B. Critical ground states in one dimension are known 
to have entanglement entropy that diverges with subsys- 
tem size for several types of pure 0, 0, M Q an d disor- 
dered d, 0, quantum critical points, with a coefficient 
determined by the central charge in the pure case. Hence, 
in one dimension, entanglement entropy gives a definition 
of universal critical entropy for quantum critical points 
that is consistent with the conventional definition (cen- 
tral charge) in the conformally invariant case. 

In higher dimensions, results have been obtained re- 
cently for gapless free fermions d, 0, and a generic 
scaling form is conjectured in A standard expecta- 
tion |lOj, which is known to be violated for free fermions, 
is of an "area law" : entanglement entropy scales as the 
area of the boundary between the subsystems. For ex- 
ample, if only one subsystem is finite, the area law pre- 
dicts that entanglement entropy scales as U 1 ^ 1 in d spa- 
tial dimensions, where L is the linear size of the finite 
subsystem. Entanglement entropy in gapped phases in 
two dimensions satisfies an area law but has subleading 
terms that probe topological order [ill ll2T] : adding and 
subtracting properly chosen regions cancels the leading 
area term and leaves a term reflecting the ground-state 
degeneracy on topologically nontrivial manifolds. 

An additional motivation for recent studies of entan- 
glement entropy in many-body physics is as a route to 
better numerical algorithms for finding ground states: 
knowing that ground states of local Hamiltonians have 



much less entanglement, even at criticality, then generic 
quantum states both explains the success of the "density- 
matrix renormalization group" algorithm in one dimen- 
sion [l?3 . IT3 . Hfl and motivates recent proposals for anal- 
ogous numerical methods in higher dimensions |16| . 

This paper obtains the entanglement entropy for the 
class of "conformal" two-dimensional quantum critical 
points 17] that includes such standard examples as the 
quantum dimer model |18l Il9| , the quantum eight-vertex 
model 17]. These quantum critical points were first in- 
troduced because perturbing away from these solvable 
critical points can yield topologically ordered phases in 
lattice models |2jj. The same method can be used to 
obtain the entanglement entropy for related noncritical 
wavef unctions. These critical points have dynamic criti- 
cal exponent z = 2 and equal-time correlation functions 
given by a two-dimensional conformal field theory (CFT). 
For the entanglement entropy created by partitioning of 
a 2D manifold into regions A and B, we find a univer- 
sal subleading logarithmic contribution in addition to the 
"area law" nonuniversal contribution proportional to the 
linear size L of the boundary: 



S = 2f s (L/a) + ac\og(L/a) + O(l). 



(1) 



Here c is the central charge of the CFT, a is the ultravio- 
let cutoff, and the coefficient a is determined by geomet- 
ric properties of the partition. The area law coefficient 
f s is interpreted below as a boundary free energy. As an 
example of the geometric dependence of a, if A is a rect- 
angle surrounded by B, a = — |, while if A has a smooth 
boundary surrounded by B, a — 0. This contribution for 
the case of a free scalar field originates in the logarith- 
mic contribution that appears in the classical problem of 
the determinant of the Laplacian operator on a 2D man- 
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ifold [2lL ( "hearing the shape of a drum" ) , which was 
extended to a general CFT by Cardy and Peschel [23|. 

The entanglement entropy of a pure state of a bipartite 
system is defined as the von Neumann entropy of the 
reduced density matrix for either subsystem: 

S = — Tr pa log pA = — Tr pb log Pb ■ (2) 

For conformal quantum critical points, the Hilbert space 
has an orthonormal basis of states | {</)}) indexed by clas- 
sical configurations {</>}, and the ground state \tpo) of the 
bipartite system is determined by a CFT action S: 

|V.o>-^L/W)e- S(W})/2 |W). (3) 
I 



Here Z c is the partition function J (d(f>) e S (W' J and ex- 
pectation values in this state reproduce CFT correlators. 

The next step is to obtain traces of powers of the re- 
duced density matrix in order to find the entanglement 
entropy using the replica trick 

d 

- Tr p log p = - lim — Tr p n . (4) 
n->i an 

First consider the reduced density matrix pa- The el- 
ement of the reduced density matrix between internal 
configurations K^i 1 }) and |{</>^}) is i after introducing a 
normalization factor to ensure Tr pa = 1, 



(W}\pa\{^}} = Tr 0B (({^}|0({^}|^)(^o(|{^})®|{^}) 

= _L I ' ^By-(S A {4>t)/2+S A {^)/2+S (<t,t,4> B )/2+S a (^,4> B )/2+S B {4> B )) _ (5) 



Here the action has been divided into regions A, B, and 
the boundary <9, where the last takes into account con- 
tributions mixing the A and B degrees of freedom (e.g., 
couplings of spins across the boundary in a lattice model) . 

Higher powers of the density matrix need not trace to 
unity: Tr p\ is now a sum over n configurations defined in 
A and n configurations defined in B. The key is to keep 
track of how these different configurations are stitched 
together at the boundary by the terms S d that link A 
and B: is linked to {<j)f } as well as {<pf +1 \ for 

i = 1, . . . , n — 1, and {4>n} 1S linked to {4>n} an( i {4>i }• 
This is normalized through division by (Z c ) n , which can 
be thought of again as n copies of A and B configurations, 
but with {4>f } linked only to {(f>f }. 

In the continuum limit, the boundary terms impose 
continuity of the fields in the CFT because strong local 
fluctuations are penalized in the action. Since each A 
field is linked to the B field of the same index, we can 
define global configurations {4>i} in both the numerator 
and denominator, but in the numerator, the links be- 
tween configuration i and configuration i ± 1 require that 
all n configurations agree on the boundary, while this re- 
quirement is absent in the denominator. Schematically, 



n Z(n configurations agreeing on the boundary) 
A Z(n independent configurations) 

(6) 

Note that this is symmetric with respect to A and B: 
Tr p^ — Tr p r g as known from Schmidt decomposition. 

If Tr p\ is known for all integer n, then assuming an 
analytic continuation, the entropy is obtained as 



<9Tr 



Pa 



On 



(7) 



The above expression for Tr p A can be put in a form that 
simplifies taking this derivative: for an explicit realiza- 
tion, consider the case of a free scalar field. Then the 
condition that n scalar fields 4>i agree with each other on 
the boundary can be satisfied by forming n— 1 linear com- 
binations ^7j(0i — ) ; which vanish at the boundary 
(i.e., satisfy Dirichlet boundary conditions), plus one lin- 
ear combination ^ J2i=i n & tnat nas n0 restriction 
at the boundary (i.e., is a free field on A U B). 

More generally, for any CFT there exists a confor- 
mal boundary condition that generalizes the notion of 
the Dirichlet boundary condition in the free scalar case. 
Thus, in terms of the partition functions Ze>, for a field 
in the whole system A UB that vanishes at the boundary, 
and Zp, for a free field in the whole system, 



and therefore 



Zj d Zjp 

Zip 



1 ZD 

Zj F 



Z F 



log 



7 A 7 B 



(8) 



(9) 



In the last equality, the Dirichlet boundary condition at 
the boundary was used to split the partition function into 
contributions from A and B, each including the boundary 
with Dirichlet boundary conditions. Finally, the entan- 
glement entropy for a general conformal quantum critical 
point is just the dimensionless free energy difference in- 
duced by the partition in the associated CFT: 



S = Fa + Fb — -F 



AUB- 



(10) 
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Explicit results for the entanglement entropy can now 
be found using results on the free energy corrections in 
conformal field theory. If AuB is the plane, A and B are 
connected (which requires that one be simply connected), 
and the boundary is smooth, then the free energy result 
of Cardy and Peschel |2^| can be applied: using region A 
as an example and supposing that it is finite, 

F A = f b \A\ + f s L- ^- log L + 0(1) (11) 

where fb and f s are the bulk and surface partial densities, 
\A\ is the area, L is the perimeter, and x is the Euler 
characteristic of the manifold 

X = 2 - 2g - b. (12) 

Here g is the number of handles and b the number of 
boundaries of the manifold. As an example, x — 1 f° r a 
planar simply connected region with boundary, such as 
a disk. In (| 1 1 p . the cutoff is now incorporated into the 
definitions of fb and f s : note that the coefficient of the 
logarithm is cutoff- independent, as changing the cutoff 
gives only an additive constant from this term. 

In the expression for S, the bulk energy terms can- 
cel, and since A U B is the plane and has no bound- 
ary, the linear in L term is just twice the dimensionless 
boundary energy in the CFT with the same regulariza- 
tion. Note that the logarithmic term differs from the 
leading (boundary) term in that it is cutoff-independent 
and hence potentially universal. The logarithmic con- 
tribution is determined for a smooth boundary by the 
change in Euler characteristic: 

c 

S = 2f s L - -(xa +XB - XAus)logi. (13) 

D 

We note immediately from this formula that if the bound- 
ary between A and B is a smooth closed curve in the inte- 
rior of AUB, there is no logarithmic contribution because 
the total Euler characteristic is conserved. Although this 
cancellation is quite specific to a smooth closed boundary, 
it shows that the existence of logarithmic contributions 
to Fa and Fb (unless they have zero Euler characteristic) 
does not imply a logarithmic entanglement entropy. 

There are two simple mechanisms that generate a log- 
arithmic contribution. The first is that the boundary 
between A and B may intersect the boundary of AU B, 
as for example if a disk is sliced into two halves, which 
modifies the total Euler characteristic. The second is that 
the boundary between A and B may have corners. [24| 
The Cardy-Peschel formula argument given here is valid 
for boundaries that are smooth except for a finite number 
of sharp corners; fractal boundaries, for example, could 
give rise to different subleading corrections. Henceforth 
we write S\ og for the possible logarithmic term in the 
entropy. Three examples are shown in Fig. 1. 

We first treat the case of sharp corners. A sharp corner 
with interior angle 7, < 7 < 2-7T, gives a logarithmic 



contribution to the free energy [22, 23] 

AF=^(l-^)logL (14) 

According to the formula (|14l> . the contribution from a 
polygon approaches that of a smooth curve as all angles 
approach ir. However, sharper angles contribute more to 
the logarithmic correction than would be predicted by the 
Gauss-Bonnet formula. Contributions from the angles 7 
and 2-7T — 7 only cancel to linear order near 7 = tt, so that 
sharp angles on the AB boundary give a net contribution. 

As an example, summing the interior and exterior 
contributions for the four corners of a rectangle A sur- 
rounded by B gives 

5io 6 = c^-i)logi = -|logi. (15) 

Note that a negative or positive sign for the logarithmic 
part of the entropy is physically permissible as a sublead- 
ing correction to the (positive) boundary entropy. 




FIG. 1: The logarithmic contribution to entanglement en- 
tropy for three partitions of a disk. Partitions (a) and (c) do 
not change the total Euler characteristic, while partition (b) 
does. Note that it is the partition length scale, which may be 
finite even if A or B is infinite, that sets R in the logarithm. 

Another source of logarithmic contributions is a 
smooth boundary that disconnects the original system. 
As an example, consider the disk of radius L cut into 
two pieces by a diameter. The resulting two half-disks 
have two sharp angles of 5 each, formed by the intersec- 
tion of the diameter with the circumference: the resulting 
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logarithmic contribution is 

Slog = 4^(1 - 4) logi = -|logL. (16) 

The quantum critical points studied in this paper can, 
as in the quantum dimer model and its generalizations 
(both Abelian and non- Abelian)|lll2Sl23, lie next to 
topologically ordered phases with subleading 0(1) terms 
in their entanglement entropy in addition to an area 
law: [ill for some nonuniversal constant a and a uni- 
versal topological number 7, 

Stopo =a£ + 7 + 0(l/£). (17) 

Another way to obtain a universal logarithmic term 
at criticality is by means of a physical process by which 
the system is split into two disjoint regions each with a 
smooth boundary. In this case there is a net change in the 
Euler characteristic leading to a lowering of the univer- 
sal term in the entanglement entropy (due to a physical 
loss of correlations) by AS — — (c/6)log£. Something 
very similar has been found to happen to the universal 
constant term 7 which arises in a topological fluid (see 
above), such as in the case of a quantum Hall fluid 27]. 

It is natural to ask whether even at criticality there are 
topological nondivergent terms in addition to the univer- 
sal logarithmic divergences studied here. For the case 
of the free scalar field, there can be nondivergent terms 
with a complicated (but scale invariant) dependence on 
the partition geometry, suggesting that there is no purely 
topological number 28]: e.g., the free energy of the annu- 
lus with radii r% < r2 contains a term loglog^/ri). 

An important question for future work is to understand 
explicitly how the entanglement evolves, for a topologi- 
cally nontrivial partition, as the Hamiltonian is tuned 
from a critical point to a massive topological phase, as 
this requires a more accurate calculation than simply 
cutting off the logarithmic divergences by a correlation 
length £ away from criticality. The dependence of the 
universal logarithmic term on the central charge c of 
the associated 2D CFT indicates that upon perturba- 
tion there should be a downward flow of the coefficient 
of this universal term. 

Although the area law at leading order is consistent 
with the general scaling ansatz of Calabrese and Cardy[3] 
in d > 1, our results show that in d = 2 there are uni- 
versal logarithmic contributions to the entanglement en- 
tropy at a class of critical points for which the "area 
law" applies. These universal contributions are deter- 
mined by the critical bulk theory and the geometric de- 
tails of the partition into subsystems. Combined with 
recent results on subleading corrections to entanglement 
entropy in massive topological phases [ill Il2j ] , it now ap- 
pears that entanglement entropy in d > 2 is considerably 
richer than the area law might suggest, even when the 
area law is applicable. 
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